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We consider a tight-binding model on the honeycomb lattice in a magnetic field. For special 
values of the hopping integrals, the dispersion relation is linear in one direction and quadratic in 
the other. We find that, in this case, the energy of the Landau levels varies with the field B as 
€n{B) rsj [(n+7)B] 2//3 . This result is obtained from the low- field study of the tight-binding spectrum 
on the honeycomb lattice in a magnetic field (Hofstadter spectrum) as well as from a calculation 
in the continuum approximation at low field. The latter links the new spectrum to the one of a 
modified quartic oscillator. The obtained value 7 = 1/2 is found to result from the cancellation of 
a Berry phase. 
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Introduction - The recent discovery of graphene has 
boosted the study of the physical properties of the hon- 
eycomb lattice, especially in a magnetic field [1]. Among 
the peculiarities of the electronic dispersion relation, 
the spectrum in the band center (e = 0) is linear, ex- 
hibiting the so-called Dirac spectrum around two spe- 
cial points at the corners K and K' of the Brillouin 
zone. Near these points the density of states varies lin- 
early [2]. In a magnetic field, it has been shown [3] that 
the energy levels around e = vary with the field B as 
e n (B) ~ ±[nB] 1 / 2 , with a degeneracy which is twice the 
degeneracy of usual Landau levels (LLs). This is easily 
explained by the two-fold valley degeneracy correspond- 
ing to the two points K and K f . This spectrum has 
to be contrasted with the familiar field dependence of 
Landau levels e n (B) = (n + \/2)eB/m for electrons in a 
quadratic band with mass m. This square root depen- 
dence has been observed experimentally [4, 5]. 

Here we present an example where the field depen- 
dence of the LLs is neither linear nor a square root, but 
reveals a new power-law, namely a [(n + j)B} 2 ^ 3 behav- 
ior, with 7 = 1/2. This is obtained for tight-binding 
electrons on the honeycomb lattice, the same problem as 
for graphene, but with special values of the hopping inte- 
grals between nearest neighbors which, contrarily to the 
case of graphene, are not taken to be equal. We find that 
around a special point of the reciprocal space, the zero- 
field spectrum is linear in one direction and quadratic 
in the other. This "hybrid" spectrum leads to a new 
magnetic field dependence of the Landau levels, properly 
described by a quartic oscillator V(X) = X 4 + 2X. The 
value 7 = 1/2 is found to result from the cancellation of 
a Berry phase. 

The model - We consider the tight-binding model on 
the honeycomb lattice with the possibility that one of the 
three hopping elements between nearest neighbors may 
take a different value t' from the two others t (fig. 1). 
This problem has been studied recently both in zero field 
[6] and in a magnetic field [7] , where the authors mention 



the special interest of the case t' = 2t. For these param- 
eters, the authors mention a square-root energy depen- 
dence of the density of states at the band center [6] , and 
they compute the evolution of the Hofstadter spectrum 
when t' varies between t and 2t [7]. Here, we emphasize 
the new and peculiar character of the LLs around e = 0, 
which has not been foreseen in previous works. 
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FIG. 1: a) Honeycomb lattice with hopping integrals t and t' , 
and elementary vectors a\ and a 2 discussed in the text. b,c) 
Low energy spectrum for the cases t' = t (b) and t' = 2t (c). 



The tight-binding Hamiltonian couples sites of differ- 
ent sublattices named A and B. The eigenvectors are 
Bloch waves of the form 
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where \Rf), \Rf) are atomic states. The sum runs over 
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vectors of the Bravais lattice. The eigenequations read 
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where a\ = a(|, a 2 



- a(|,— -^) are elementary 
vectors of the Bravais lattice, a is the interatomic dis- 
tance, and t, are shown in fig. l.a. It is well known that 
when t f = t, the energy vanishes at two points D and D' 
which are located at the corners K and K' of the Bril- 
louin zone (K = 2a*/3+aJ/3, K 1 = e^/3+2a5/3, where 
a* and a\ are reciprocal lattice vectors). As t' increases, 
the two points D and D' approach each other (their dis- 
tance varies as DD' = -KK' arctan \J ^t 2 /t' 2 — 1) and 
merge into a single point Do = (a* + a^)/2 when t' = 2t 
(for £' > 2t, a gap opens between the two subbands). An 
expansion k = D + p around these two points gives the 
low-energy spectrum, solution of the two equations (here 
a = 1 for shorter notations) 



tck = [-lit'p x ±V3t"p y + ^t'(3pl+p 2 y ) + ^-t"p xPy ]4 



with t" = \/t 2 — t' 2 /4 and a similar expression for ec^ 
with (z — > — i). The ± sign denotes the vicinity of the 
two points D and D' . 

We now consider the case t r = 2t. In the vicinity of 
the single point Do, keeping the leading terms in each 
direction, we find 
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eci = ta(-Sip x + -ap 2 y )c% 

and a similar equation with (A —> B,i — >• - 
the Hamiltonian can be written in the form 
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-i), so that 
(4) 



where we have defined the velocity c = 3ta, an effective 
mass m = 2/(3ta 2 ), and a "mass energy" mc 2 = 6t (we 
fix h = 1). The eigenvalues read 
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Remarkably, the spectrum is linear in one direction, 
quadratic in the other (fig. l.b,c). It has been noticed 
[6] that such a dispersion relation leads to a square root 
dependence of the density of states, that can be written 
in the form p(e) oc ^^^/e which is unusual in two dimen- 
sions. One may expect that this peculiar behavior leads 
to a new repartition of LLs. 

Effect of the magnetic field, new Landau levels - To 
describe qualitatively the effect of a weak magnetic field, 



we first start with a simple semiclassical argument. The 
quantization condition for energy levels in a magnetic 
field B has the form S(e) = 27r(n + j)eB where 5(e) is 
the area of an orbit of energy e in reciprocal space and 
7 is a constant < 7 < 1 [8]. From eq. (5), one finds 
easily [9] 



c 



\B} 2 ^ 
(6) 

with (3 = 2r(l/4) 2 /(3v / 7r) ~ 4.9442. This new behavior 
has to be contrasted with the usual case of free massive 
particles where 5(e) = 27rrae, so that e n = uo c (n + 7) 
or with the case of Dirac particles where 5(e) = 7re 2 /c 2 , 
leading to a square root magnetic field dependence e n = 
±c^2eB(n + 7) of the energy levels. The phase factor 7 
cannot be obtained from such semiclassical argument. 

The B 2 / 3 dependence can also be obtained from an- 
other simple argument related to the new energy depen- 
dence of the density of states. Since it has a square root 
dependence near e = 0, one may expect a new field de- 
pendence of the LLs, simply because the repartition of 
the zero field states into LLs is different. Assume quite 
generally a power-law dependence of the zero field den- 
sity of states, p(e) ~ e K and a given number n of filled 
LLs. The total number of states in these levels, neB/h, 
has to be equal to the total number of filled states which 
is or order p(e)de so that e n (B) ~ [(n + 7)5] «+i 

[10]. A constant density of states leads to the usual LLs, 
a linear density of states to the square-root behavior ob- 
served in graphene. Here, the same argument gives a new 
power law e n (B) ~ [(n + j)B]s for the LLs. 

Fig. 2 presents the spectrum of the tight-binding prob- 
lem in a magnetic field (the so-called Hofstadter spec- 
trum) for the honeycomb lattice. It has first been calcu- 
lated by R. Rammal for the case t' — t and more recently 
for t' = 2t [7, 12]. The procedure to obtain this spectrum 
is described at length in ref. [12]. The fractal structure 
results from the competition between magnetic field and 
lattice effects. For a commensurate reduced flux / = p/q, 
the spectrum exhibits 2q sub-bands. It is a periodic func- 
tion of the reduced flux / = Ba 2 3^/(200) through one 
plaquette in units of the flux quantum 0q. At low field, 
the lattice effects are negligible, and we expect to recover 
the results of a continuum limit. The linear dependence 
of the LLs is clearly seen on the top and bottom of the 
spectrum (fig. 2). For t' = t, the square-root dependence 
of the levels is observed around e = 0. For t' = 2t, the 
levels exhibit the new magnetic field dependence which is 
well fitted by a power law B 2 / 3 (fig. 2.d). In order to de- 
rive the low-field spectrum analytically, we now calculate 
the magnetic field effect in the continuum approximation 
around e = 0. 

Spectrum in the continuum approximation - We now 
apply a magnetic field B and use the Landau jauge A = 
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FIG. 2: a) Hofstadter-Rammal spectrum in the case t' — t. 
b) Spectrum for the case t! = 2t. c) Low-field behavior in the 
case t! — t. The low-field bands are well fitted by the analytical 
expansion e n = ±t\/27r3 1 / 2 nf . d) In the case t' — 2t, the low- 
field spectrum is well fitted by eq. (12). Deviations at higher 
field or higher energy are due to lattice effects. 



(0, Bx, 0). The substitution p y — > p y — eBx leads to the 
new Hamiltonian 
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with an effective cyclotron frequency uo c = eB/m = 
3eBta 2 /2 which can be also written in terms of the re- 
duced flux / = 5a 2 3v / 3/(20 o ) through one plaquette of 
the lattice, namely uj c = ^tf. In x = x — p y /eB, quan- 
tization of p y leads to the usual degeneracy of Landau 
levels. The energy levels e are solutions of 



(c 2 P 2 x + (\™u 2 c ) 2 x A - i^m^ 2 c [p x ,x 2 ^ = e 2 ^ . 

Introducing new dimensionless conjugate variables X and 
P, we rewrite 

/ 2 2 \ 2 / 3 

(^f-j (P 2 +X 4 -z[P,X 2 ])^ = e V. (7) 



This expression shows that the eigenvalues necessarily 
scale as B 2 ! 3 . Neglecting first the linear term [P, X 2 ] = 
—2iX, the eigenvalues of the quartic oscillator P 2 + X 4 
can be easily estimated, at least for large n, in the WKB 
approximation [13] and are found to be of the form 
C(n + 1/2) 4 / 3 with C = [3^V^/r(l/4) 2 ] 4 / 3 ~ 2.18507. 
Therefore, the Landau levels are given by (restoring K) 



±A(mc 2 ) 1/3 l(n + l/2)hu; c } 2/3 



(8) 



with A = y ^3 — 1-17325. This is precisely the depen- 
dence expected from the above semiclassical argument 
(6), with the phase factor 7 now determined to be 1/2. 
Replacing m and uo c by their expressions in terms of the 
lattice parameters considered here, we finally obtain 



±a[(n + l/2)/] 2 / 3 



(9) 



where a = (2tt) 2 / 3 v / C ~ 5.0333. The f 2 ' 3 dependence 
of these levels is clearly seen in fig. 2.d, and their (n + 
1/2) 2 / 3 dependence is confirmed on fig. 3. 




FIG. 3: For / = 1/g, plot of qe^ 2 vs n (*), and comparison 
with the WKB solution (9) of the quartic oscillator (straight 
line) and with the numerical solution of the modified quartic, 
X A + 2X oscillator (o) for n = 0, 1. For n > 2, the two 
solutions are indistinguishable. 

It is interesting to compare eq. (7) with the eigenequa- 
tions for Dirac fermions: 



eBc 2 (p 2 + X 2 -i[P : X]^ = e 2 ip , 



(10) 



with e n = ±cy2ePn or for free massive particles 

(P 2 +X 2 )^ = 6^, (11) 
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with e n = (n+ l/2)ftw c , and to emphasize the role of the 
commutator which enters eqs (7,10). There are several 
discussions in the literature [14] to explain the disappear- 
ance of the phase term 1/2 in the case of the Dirac spec- 
trum, and we return to this point later. We notice from 



eq. (10) that this disappearance comes simply from the 
commutator [P,X] = —i. In our case, the commutator 
which appears in eq. (7), [P, X 2 ] = — 2zX, modifies the 
quartic potential which becomes X 4 + 2X. This linear 
term is actually a small perturbation negligible when n 
is large. Taking into account this linear term, a numeri- 
cal calculation of the eigenvalues of this modified quartic 
oscillator finally gives (fig. 3) 

e n = ±ag{n)[[n + 1/2) ff 3 (12) 

where g(0) ~ 0.808,0(1) ~ 0.994. For n > 2, g(ri) ~ 1 
so that the WKB solution (9) of the quartic oscillator 
turns out to be extremely good. Finally, let us mention 
that since the two Dirac points D and D' have merged 
into a single point Dq for t' = 2t, the valley degeneracy 
has disappeared and the LLs degeneracy has recovered 
its usual value. 

We now comment on the relation between the phase 
factor 7 entering the quantization of semiclassical orbits 
and a Berry phase. It has been established that [14] 




where V is the contour of a semiclassical orbit. It is easily 
found from eqs. (1,2) that, for the isotropic case t f = t, 
§ v (k\Vh\k) 'dk = ±7r around the points D and D' which 
can be seen as topological defects (see fig. 4. a), leading to 
7 = as well known for graphene. For our case t' = 2t, 
these two topological defects D and D merge into a 
single one Do and annihilate. As may be seen on fig. 
4.b, we now obtain f r (k\Vk\k) • dk = 0, which explains 

TT rl^l T T ^ — 1 /O 



sion relation. Recent experiments and calculations for 
graphene have shown that a linear dispersion relation 
leads to a square-root dependence of the LLs with the 
field. In this paper, we have presented an example with 
an usual dispersion relation, which leads to a square-root 
energy dependence of the density of states and to a new 
field dependence \{n + 1/2)P] 2/3 of the LLs. More gen- 
erally, we may consider a dispersion relation of the form 
e = (p% +Py) 5 . The area of the orbits of a given energy e 

a + P 

varies as S(e) oc e <*p s , so that the Onsager quantization 
rule for energy levels leads to the general dependence of 
the Landau levels with the magnetic field 

a/35 

e n (B) ~ [(n + 7)£]^F* . 

This dependence can also be obtained from a counting 
argument based on the energy dependence of the density 
of states which is easily found to vary as p(e) ~ e^ + ^ _1 . 

To conclude, we briefly comment on possible exper- 
imental realizations leading to such an electronic spec- 
trum. An anisotropic version of the graphene structure 
with t' 7^ t, called quinoid and discussed long ago by 
L. Pauling [15], could be induced by uniaxial stress or 
bending of a graphene sheet. We also mention recent dis- 
cussions on the feasibility of such a structure with cold 
atoms in an optical lattice created by laser beams [16]. 
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FIG. 4: Wave-vector dependence of the ratio cf /c£ for the 
eigenvectors (1). The figures represent a part of the Brillouin 
zone near the points K and K' . a) For t' — t, the Berry 
phase ±7r on semiclassical orbits around the points D — K , 
D' — K' implies 7 = 0. For t! — 2t, the Berry phase on a 
contour around Do vanishes, implying 7 = 1/2. 

Conclusion - The field dependence of the LLs depends 
dramatically on the structure of the zero field disper- 
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